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Abstract. We study theoretically the phoretic motion of a spheroidal particle, which generates solute 
gradients in the surrounding unbounded solvent via chemical reactions active on its surface in a cap-like 
region centered at one of the poles of the particle. We derive, within the constraints of the mapping to 
classical diffusio-phoresis, an analytical expression for the phoretic velocity of such an object. This allows 
us to analyze in detail the dependence of the velocity on the aspect ratio of the polar and the equatorial 
diameters of the particle and on the fraction of the particle surface contributing to the chemical reaction. 
The particular cases of a sphere and of an approximation for a needle-like particle, which are the most 
common shapes employed in experimental realizations of such self-propelled objects, are obtained from 
the general solution in the limits that the aspect ratio approaches one or becomes very large, respectively. 
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1 Introduction 

The increasing interest in the development of "lab on 
a chip" devices and of drug-delivery systems has led to 
a stringent need of scaling standard machinery down to 
micro- and nano-scales. This reduction in length scale has 
raised a number of challenging issues, such as developing 
ways to enable small objects to perform autonomous, di- 
rectional motion [T][2]- 

Although the experimental and theoretical research in 
this area is still in its early stages, several proposals for 
such "self-propellers" have already been tested experimen- 
tally (see, e.g., Rcfs. |Il[2j[3l[4l[5]); a review of the recent 
progress in this field can be found in Ref. [Hj . These "proof 
of principle" proposals have generally employed particles 
with axial symmetry, i.e., cylindrical rods [TJ[5] or spheres 
[HIS], because they are relatively easy to manufacture, al- 
low for a good control of the desired surface modifica- 
tions, and their simple geometry is a significant bonus 
for the theoretical analysis of the experimental results. 
The underlying idea, as put forward by Whitesides and 
co-workers [T], is that an asymmetric decoration of the 
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particle with a catalyst, which promotes an activated reac- 
tion in the surrounding liquid medium generating product 
molecules, can provide motility through a variety of mech- 
anisms. As in the initial design proposed in Ref. [T], the 
simplest example is the propulsion of mm-size objects due 
to the ejection and subsequent bursting of bubbles formed 
by the product molecules [O2 for PDMS plates with Pt 
catalyst tails placed in hydrogen peroxide (H2O2) aque- 
ous solutions]. As the size of the particle is decreased to- 
wards the micron scale or below, viscous and surface forces 
start to dominate and inertia-based mechanisms such as 
the "bubble ejection" propulsion become ineffective. If the 
product molecules remain dissolved in the surrounding liq- 
uid medium as, e.g., in the experiments reported in Refs. 
[HI11[S| , the result of an asymmetric distribution of catalyst 
is that the chemical reaction gives rise to concentration 
gradients along the surface of the particle. It has there- 
fore been argued [51[7] that in such cases the motion of the 
catalyst-covered "active" particle is rather phoretic, i.e., 
the result of the interactions between the particle and the 
non-uniformly distributed product molecules generated by 
the chemical reaction. For example, in the case of Au- 
Pt rods in H2O2-H2O mixtures [5] the product molecules 
O2 play the role of a solute the concentration gradient of 
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which in the solution (formed by the H2O2 - H2O mixture 
as the solvent and O2 as the solute) is the field that may 
induce phoretic motion. 

In many cases the magnitude and the direction of the 
experimentally observed phoretic velocity of such "active" 
particles are compatible with a variety of microscopic mech- 
anisms, such as surface tension gradients [UinilH] (note that 
Ref. [3| provides an elegant example that rotational mo- 
tion can also be achieved), cyclic adsorption and desorp- 
tion [5], electrokinetics [Hini!l!0, or diffusio-phoresis [31[71 
[TOlIll) . (By using the notion of diffusio-phoresis, here we 
refer strictly to phoresis due to gradients of a solute, i.e., 
we do not consider the case of bimetallic particles [HE] for 
which charge transfer and electro-chemistry may be the 
dominant effects. For simplicity, we focus here on the case 
of electrically neutral solutes, which in the literature is 
often also called "chemo-phoresis".) Thus understanding 
these systems and discriminating between these various 
possibilities require a careful, detailed theoretical analysis 
to predict the dependence of the velocity on the control 
parameters of the system, such as, e.g., the H2O2 con- 
tent of the aqueous hydrogen peroxide solvent used in the 
experiments reported in Refs. [USKS] or the fraction of 
the particle surface which is catalytically active. With the 
notable exception of Ref. [10], which has used a micro- 
scopic description for the interfacial region at the expense 
of having to carry out numerically most of the analysis, 
the theoretical approach so far has been to map these 
systems onto the case of classic phoresis in an externally 
imposed gradient of a field such as, e.g., a solute concen- 
tration or an electric potential (see Ref. [12] for a review of 
the theory of classic phoresis and additional references; a 
detailed discussion of the shortcomings of such a mapping 
is provided in Ref. ^^). 

Starting from the model system proposed in Refs. [H 
[7] , here we study the phoretic motion of a spheroidal par- 
ticle which generates number density gradients of product 
molecules emerging from chemical reactions, which are ac- 
tive on the surface of the particle in a cap-like region cen- 
tered at one of the poles. The product molecules diffuse 
into the surrounding unbounded three-dimensional New- 
tonian liquid solvent. Similar to the earlier studies in Refs. 
|71 [TT1[T^[H] ■ our work is based on adopting the standard 
theory of phoresis for the present case, in which the gra- 
dients are self-generated rather than being produced and 
maintained by external sources |15j . The motivation for 
this work is to provide a unified description (within the 
standard theory) for the diffusio-phoretic motion of ob- 
jects belonging to an extended class of geometrical shapes 
relevant to experimental studies [^l^[51[TB] . The previously 
studied spherical and needle-like shaped objects are re- 
covered as particular limiting cases. We note that here we 
focus on the case of rigid particles. If the body is actually 
soft and deformable various additional phenomena, such 
as a transfer between translational and rotational motion 
upon shape changes, may occur (see, e.g., Ref. [IT])- 

The outline of the paper is as follows. In Section [2] we 
define the model. Section [3] is devoted to the derivation of 
the diffusio-phoretic velocity; it includes also the compu- 



tation of the distribution of the product molecules which 
induces the phoretic motion. The results, as well as the 
connections with the previous studies in Refs. I^ITlfTTlIB]. 
are discussed in SecH] We conclude with a brief summary 
in Section [5] 



2 The Model 

The system we consider is shown in Fig. [TJa). It consists 
of an impermeable, spheroidal, rigid particle of polar and 
equatorial semi-axes i?i and R2, respectively. At one of the 
poles there is a cap-like region (the black area in Fig. [IJ 
covered by a catalyst [with density a (number of catalyti- 
cally active sites/unit area)] which promotes the chemical 
conversion of a surrounding solvent (not shown in Fig. 1) 
into product molecules of diameter a [small hatched cir- 
cles in Fig.[IJa)]. (This is a so-called "Janus particle" [TT].) 





-0.87 r 



Fig. 1. (a) An impermeable, spheroidal particle of polar and 
equatorial semi-axes Ri and R2, respectively, with a cap-like 
part of the surface covered by a catalyst (depicted as a black 
area). The aspect ratio of the particle is Sr = R2/R1', Sr < 1 
(shown here) corresponds to a prolate shape, while Sr > 1 
to an oblate one. Sr — 1 corresponds to a sphere. The prod- 
uct molecules of diameter a are shown as small hatched cir- 
cles. O denotes the geometric center of the particle from which 
the number density p(r, t) of product molecules is measured, 
(b) and (c) show cuts of the xz plane through the prolate 
and oblate iso-surfaces in terms of prolate (^, 0, 7;) and oblate 
(CfptX) spheroidal coordinates, respectively [see, c.f., Eqs. 
(fT2)l and Eqs. (f20ll - ((24)l . respectively]. ^0 and (0 denote the val- 
ues of the prolate ^ and oblate coordinates, respectively, for 
which the corresponding iso-surfaces coincide with the surface 
of the particle (shown as a thick gray line; the thick black line 
at the lower pole indicates the catalyst covered region). 



The shape of the object is characterized by its aspect 
ratio Sr = R2/Ri'. Sr < 1 refers to a prolate spheroid, 
Sr > 1 refers to an oblate spheroid, and Sr = 1 to a sphere 
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Fig. 2. A schematic representation of the spheroidal shapes 
and of the limiting cases which are discussed in the main text; 
Sr ~ R2/R1 is varied by varying R2 while keeping _Ri fixed. 

(see Fig. [5]). The various spheroidal shapes, i.e., the val- 
ues of Sr in the full range < < 00, are systematically 
explored by considering the parameter Ri fixed (introduc- 
ing a characteristic length scale) while varying the param- 
eter i?2- The cartesian coordinate system is chosen such 
that for prolates the foci of the generating (through ro- 
tation around its major axis) ellipse are located on the 
z axis symmetrically with respect to the origin O which 
is located at the geometrical center of the particle; corre- 
spondingly, for oblates the foci of the generating (through 
rotation around its minor axis) ellipse are located in the 
xy plane symmetrically with respect to the origin O. The 
ratio Sh = h/Ri between the height of the cap- like cata- 
lyst covered area and the polar semi-axis characterizes the 
fraction of the particle surface covered by catalyst. Note 
that the case of a particle with a point-like catalytic site 
located at one of its poles can be included in the model 
by considering the limit {sh — >■ 0, a — >■ 00} under the con- 
straint that cr X Sh is finite (see, c.f., Appendix IB|) . There- 
fore, by varying the parameters Sr and Sh the geometries 
studied in the literature are recovered as the following lim- 
iting cases (see also Ref. [IHl (a)]): 

— {sr ^ 1~,0 < Sh < 1} (where 1~ indicates that the 
limit is taken through prolate shapes, i.e., Sr ^ 1) cor- 
responds to a sphere partially covered by catalyst as 
in Ref. [TT]. (Note that it can be similarly obtained 



through oblate shapes, i.e., ^ 1, as the limiting case 

{Sr ^ l+,0 < Sft < 1}); 

— {sr — > l^,Sh — >■ 0} (or, equivalently, {sr — > ^ 
0}) with ash finite corresponds to a sphere (obtained 
as the limit of prolate or oblate shapes, respectively) 
with a point-like catalytic site as in Ref. [7]; 

— {sr ^1, < s/i < 1} corresponds to the limit of a 
needle-like particle, which approximates an elongated 
cylinder (rod) partially covered by catalyst as consid- 
ered in Refs. WEEl- 

Additionally, 

^ {sr 3> 1, < Sh < 1} approximates a thin disk par- 
tially covered by catalyst. 

In general, the chemical conversion of a solvent gen- 
erates several types of product molecules. Here we shall 
focus on the particular case in which the chemical conver- 
sion A'^ A' + B of a, solvent molecule {A) leads to two 
molecules {A' and B) only, one very similar in size and 
properties with the solvent itself {A' A), the other one 
{B) significantly different. In the following only this latter 
is denoted as "product molecule" and plays the role of a 
solute in the solvent. In other words, we consider a situa- 
tion in which the net result of the chemical conversion can 
be approximated as the generation of a solute, solely, and 
in which the reaction does not lead to a solvent depletion 
near the catalytic site (which otherwise would act as a 
solvent sink). For example, this is approximately the case 
for the Ft catalyzed decomposition of hydrogen peroxide 
(H2O2) in aqueous solution into water (II2O) and oxygen 
(O2) molecules, as discussed in Refs. 013]. In these exper- 
imental studies the oxygen plays the role of the product 
molecule the properties of which differ significantly from 
those of the solvent. Here the solvent is actually a binary 
liquid mixture of H2O and H2O2 for which H2O is chem- 
ically passive and does not participate in the chemical 
conversion. 

We thus assume that the reaction at the catalytic zone, 
i.e., the cap-like area centered around the pole at —Riez 
(where is the unit vector of the z-axis), acts effec- 
tively only as an ensemble of independent sources - uni- 
formly distributed over the cap area [with number den- 
sity fT (number of catalytically active sites/unit area)] - 
of product molecules of diameter a, which are diffusing 
into the solvent with diffusion coefficient D [IT]. We shall 
focus on the case in which the reaction rate vb at a cat- 
alytic site, i.e., the number of product molecules created 
per unit time, is independent of time. The number den- 
sity of product molecules is considered to be so low that 
among themselves they behave like an ideal gas. There 
is an interaction potential between the product molecules 
and the moving particle, which includes the impermeabil- 
ity condition at the surface of particle. The interactions 
between the product molecules and the solvent are ac- 
counted for in an effective way via the Stokes - Einstein 
expression D = ksT / [Ziriia) for the diffusion coefficient 
D of the product molecules [12], where ks is the Boltz- 
mann constant, T is the temperature, and /x is the viscos- 
ity of the solution (solvent plus the solute, i.e., the product 
molecules) . 
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3 Phoretic velocity 

3.1 General considerations 

The presence of a source of solute (product molecules) on 
parts of the surface of the particle creates a non-uniform 
and time dependent distribution of solute in the solution 
(see Fig. [1}, i.e., a non-uniform composition of the solu- 
tion. The solute number density p(r, t) is characterized by 
two important features. On the scale of the particle size 
min(i?i,i?2) ^ a, p varies due to the diffusion process. 
On the much smaller length scale of the solute diameter 
a, p{r, t) varies also near the particle surface because it 
interacts with the particle via an effective substrate po- 
tential ^ (in the sense that W describes the interaction of 
the solute molecule with the particle in the presence of 
the solvent). Typically the range 5 oi^ \s proportional to 
the solute diameter a. Accordingly the solute molecules 
interact directly with the particle only if they are within 
a thin surface film of thickness (5, which is assumed to be 
not deformed by the motion of the particle [12 . 

In what follows, we assume that the solvent can be 
considered as a continuum with constant density psoiv 
even at the length scale 5, such that a hydrodynamic de- 
scription applies within the aforementioned surface film. 
Within this picture the solute molecules and their effective 
interaction with the particle are replaced by a correspond- 
ing distribution of "point forces" acting on the solvent 
in the film. Within the limitations of such an approach 
the hydrodynamic description of the solution naturally 
splits into that of an "inner" region formed by the surface 
film and that of an "outer" region formed by the exterior 
space beyond the surface film. Furthermore we assume 
that for typical velocities V in phoresis both the Reynolds 
number Re ~ psoivV Taa,yi{Ri, R2) / p., where Psoiv is the 
mass density of the solvent, and the Peclet number Pe ~ 
y max(i?i,i?2)/-D are small, such that one can approxi- 
mate the hydrodynamic description with the creeping flow 
(Stokes) equations and disregard the convection of the so- 
lute compared to its diffusive transport. (Here we have 
assumed that the magnitude of the hydrodynamic flow 
velocity u is similar to that of the phoretic velocity V , an 
assumption which will be justified a posteriori.) For a par- 
ticle of size max(i?i,i?2) = 10 pui moving through water 
(density p = 10^ kg/m'^, viscosity p = 10""^ Pa s) with a 
velocity of the order of puv/ s, which is typical for phoresis, 
one has Re ~ 10~^. For the diffusion at room tempera- 
ture {ksTroovi 4 X 10-21 J) of O2 (a - 10-^° m) in 
an aqueous H2O2 solution {p ~ lO"'^ Pa s), the Stokes- 
Einstein relation leads to an estimate D ~ 4 x 10"^ m^/s 
for the diffusion coefficient (in agreement with Ref. [B]), 
and thus Pe ~ 10^^. Therefore both latter assumptions 
are justified. 

Following Refs. [T^[^[TT| . the ensuing asymmetric, 
non-uniform solute number density p{r, t) around the par- 
ticle will give rise, within the surface film only, to a pres- 
sure gradient along the surface of the particle. This is the 
case because the surface film is very thin on the scale of 
the particle size min(i?i,i?2) so that the equilibration of 
the composition profile of the solution within the surface 



film can be assumed to be fast along the direction normal 
to the surface of the particle compared with the diffusional 
relaxation time of the composition gradient along the sur- 
face of the particle, which typically involves a length scale 
of the order of max(i?i , i?2). In accordance with our earlier 
assumption that the solute particles can be considered to 
form an ideal gas, this implies that near the surface of the 
particle the spatial variations of p (within the surface film 
and in the direction normal to the surface of the particle) 
are given by a Boltzmann distribution exp(— /?!?'), where 
/? = l/(fcBT), corresponding to the local equilibrium con- 
figuration in the presence of the effective interaction po- 
tential between the particle and the solute molecules, with 
a prefactor which depends on the position along the sur- 
face of the particle [T2|. Mechanical equilibrium of the 
solvent within the surface film along the direction normal 
to the surface of the particle (no flow along this direc- 
tion) requires the pressure gradient along the normal to 
be equal to the body force densities due to the effective 
particle-solute interactions. Therefore the pressure within 
the surface film differs from the "outer" pressure field by 
an "osmotic pressure" term, i.e., a term proportional to 
the extra solute density Ap in excess to that density in the 
outer region. Since this osmotic pressure varies along the 
surface of the particle, there is a pressure gradient along 
the surface of the particle. This lateral pressure gradi- 
ent is not balanced by any body force and thus generates 
shear stress within the surface film. Therefore it induces 
hydrodynamic flow of the solution along the surface of the 
particle and entails motion of the particle with a steady- 
state velocity V pTll^ . Because the system has azimuthal 
symmetry, the motion is along the z-axis, i.e., V = Vez. 

We note that in a number of instances a different ar- 
gument has been used in deriving an expression for the 
velocity of the particle. It invokes a balance between a 
Stokes-likc viscous friction and a "driving force" produced, 
e.g., by a spatially non- uniform surface tension (owing to 
the product molecules changing the composition of the so- 
lution) [31[5]. Such an argument (see, e.g., Ref. [5] and the 
follow-up Ref. [5] which aimed at improving the descrip- 
tion through a more rigorous calculation of the density 
of product particles) is in sharp contradiction with the 
generally accepted view that in phoresis the motion of 
the particle plus the surface film (i.e., the region within 
which the interaction between the particle and the solute 
molecules is relevant) is force-free, with the velocity be- 
ing determined precisely by this requirement [T^[TT1[T51 
[H] . This confusion leads to results which differ from the 
correct ones roughly by a factor which is the square of 
the ratio between the size of the particle and the range of 
the particle-solute interactions; i.e., the predicted values 
of the velocity are too large by orders of magnitude. 

3.2 Phoretic slip and phoretic velocity 

As discussed above, the pressure gradient along the surface 
of the particle, induced by its interaction with the non- 
uniform distribution p{r, t) of product molecules within 
a thin surface film, leads to fiow of the solution relative 
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to the particle. Assuming that p(r, t) is changing slowly 
in time, this hydrodynamic flow, considered to be locally 
planar, translates into a (phoretic) slip-velocity, 

Vs(rp) = -&V^p(rp) , for Vp & Us , (1) 

as a boundary condition for the hydrodynamic flow in the 
outer region (T^UHIIISI ■ In this equation Eg denotes the 
outer edge of the surface film (which is a surface at a 
distance S parallel to the surface S of the particle), 
denotes a point P on Sg, denotes the projection of the 
gradient operator onto the corresponding local tangential 
plane of the surface of the particle (actually of Us; but for 
the outer problem this can be replaced by S because the 
variations of p are over length scales much larger then 6) , 
while 

6 = (2) 

is an effective "mobility" coefficient, and A = \/\A\ is a 
characteristic length scale. The latter is given in terms of 
the effective interaction potential iZ' between the particle 
and the product molecules. Within a local coordinate sys- 
tem in a small domain of the surface film of width S, cen- 
tered at position r^, !Z' determines the product molecule 
distribution along the direction y normal to the particle 
surface and yields [T^[T5] 

oo 

A^Aivp)^ J . (3) 



We have explicitly indicated in Eq. ^ that the effec- 
tive interaction potential \P and, consequently, the length 
scale A may vary slowly along the surface of the particle 
over length scales much larger than the thickness 6 of the 
surface film. This is the case because one may reasonably 
expect that the effective interactions between the particle 
and the product molecules depends on the local chemical 
composition of the surface. Thus they can be different in 
the region Sc at the pole which is covered by the cata- 
lyst from that at the chemically inert part of the particle 
surface. We note that such non-uniformity of the surface 
properties has been also explicitly noted in, e.g., Ref. [2], 
anticipating that the surface tension of the rod-solution 
interface is different at the Pt and Au ends of the rod. 
However, its role was ignored in the analysis there, as well 
as in the follow-up work in Ref. [8], based on the argu- 
ment that the density of O2 is uniform over the Pt end. 
Nevertheless, such an argument is clearly contradicted by 
the theoretically calculated O2 density along the surface of 
the rod which shows significant variations over the whole 
surface of the rod (see Fig. 4 in Ref. [5]). In a first or- 
der approximation one can account for these non-uniform 
properties along the surface by using a position dependent 
effective mobility b{rp) [Eq. (0) with A — >• ^(r^,)], which 
can attain two values, in the expression of the phoretic 
slip- velocity [Eq. ([T])]. This approach has been pursued in 
Ref. [11]. In the case of a sharp boundary between the 
catalyst-covered region and the inert one such an approx- 
imation is probably justified. However, if the non-uniform 



properties vary smoothly and significantly over the sur- 
face, then it becomes unclear if this can be accounted for 
solely by a piece-wise constant effective mobility b{rp). In 
such detailed analysis of the hydrodynamic flow 

in the surface film forming the inner region is required 
for understanding the effects of inhomogeneities on the 
induced phoretic motion. For example, it has been shown 
that a solid sphere with position-dependent slip boundary 
condition immersed in a laminar flow experiences a torque 
induced by this inhomogeneity |23| . Moreover, in the case 
of electrophoresis it has already been shown that for a par- 
ticle with non-uniform surface properties the phoretic mo- 
tion depends strongly on the details of this non-uniformity, 
which may give rise to counter-intuitive results such as 
electrophoresis of spherical particles which are electrically 
neutral and have a zero mean (^-potential (see, e.g., Refs. 
[2^I25| ). But there may be cases in which such coatings 
by certain catalysts do not significantly change the ef- 
fective interaction between the particle and the product 
molecules. Accordingly, here we shall focus on the case in 
which <]/, and thus the effective mobility b, can be assumed 
to be constant over the surface of the particle. The general 
case of an effective interaction potential which varies over 
the surface of the particle will be discussed elsewhere. 



As mentioned above, the dynamics at small Pc num- 
bers implies that the convection of the solute is negligible 
compared to the diffusive transport, i.e., the diffusion of 
the solute is apparently decoupled from the hydrodynamic 
flow. However, a coupling between the density profile of 
the solute, i.e., the solute diffusion, and the flow of the so- 
lution is re-established by Eq. ([T]) which, in the reference 
frame co-moving with the particle, represents the bound- 
ary condition (BC) at the edge of the surface film for hy- 
drodynamic flow in the outer region. Because there are 
no forces acting on the solution beyond the surface film, 
based on the assumption of low Re numbers the hydrody- 
namic flow in the outer region is obtained as the solution of 
force free and incompressible Stokes equations subject to 
the following BCs in the laboratory frame: (i) prescribed 
velocity V-t-v^ at the edge Us — E of the surface film (i.e., 
sticking on the surface of the particle plus a slip velocity 
Vs at the edge of the surface film) and (ii) zero velocity 
(fluid at rest) far away from the particle. After computing 
the solution in the outer region, which depends paramet- 
rically on the velocity V of the particle via the above BC 
(i), the phoretic velocity is determined by the condition 
that the motion of the system composed of the particle 
plus the surface film is force free. Here the cumbersome 
explicit calculation of the hydrodynamic flow is avoided by 
using Brenner's generalized reciprocal theorem |18[ (b)], 
which allows one to express the phoretic velocity of the 
particle as a surface integral of the phoretic slip velocity 
weighted by the normal component of the position vector 
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(see Appendix A and, e.g., Ref. }25j'): 



V 



3V„ 



(fi • rp) (e, • V^p) dS, 



(4) 



where the second equality follows by using Eq. (P) to re- 
place Vs. 

In this equation Vp = (4/3)7rs^i?f denotes the volume 
of the particle, n is the unit vector of the direction nor- 
mal to the surface of the particle, and in the integral we 
replaced Ss by S because the flow in the outer region 
varies over length scales which are much larger than the 
difference S between the two. According to Eq. (j?]), knowl- 
edge of the number density p of product molecules at the 
outer edge Us of the surface film completely determines 
the velocity V of the particle. 



3.3 Steady-state distribution of product molecules 
(solute) 

After switching on the catalytic reaction there is an ini- 
tial transient phase during which the distribution p{r, t) of 
product molecules builds up. Because we are dealing with 
an unbounded solvent, one may expect that at long times 
(t — > oo) the distribution of product molecules reaches a 
steady-state distribution p{y). In what follows, we focus on 
the motion in the long-time regime (in practice, at times t 
much longer than the diffusion time to = [max(i?i , i?2)]^/-D 
for the numerical estimates from Subsec. 13.11 to ~ 2.5 s) 
in which p and consequently the phoretic slip Vg and the 
phoretic velocity V [see Eqs. ([T]) and (HJ, respectively] be- 
come time independent. Note that here we implicitly as- 
sume that the time-scale ta of the rotational diffusion of 
the polar axis of the particle is sufficiently larger than to 
such that the steady-state distribution piv) is attained and 
translational motion of the particle with uniform velocity 
V occurs; at times t ^ tn the rotational diffusion of the 
polar axis leads to a quasi-diffusive behavior of the parti- 
cle displacement (see, e.g., Refs. [4l[26]). These two regimes 
are clearly observable in the experiments discussed in Ref. 
[2] (see, in particular, the supporting information therein). 

Within the assumptions that the diffusion of product 
molecules is fast compared with the convection by the sol- 
vent flow, i.e., in the limit of small Peclet numbers, and 
that the product distribution p(r) in the steady state is 
undisturbed by the flow, i.e., neglecting any so-called po- 
larization effects of the surface film [T^], the steady state 
distribution p(r) of product molecules in the outer region 
around the moving particle is governed, in the co-moving 
frame, by the diffusion equation 



film and (ii) at each point r Cz Ss the product molecules 
current in the outward direction n normal to the sur- 
face is equal to the total reaction rate in an infinitesimal 
element of S centered at r — i5n. The latter BC holds 
within the assumption that the surface film is very thin 
[S ^ min(i?i, i?2)] such that in the steady state the lat- 
eral transport of product molecules in the surface film is 
negligible compared to the transport into the outer region 
along the direction normal to Ss- Since in the calcula- 
tion of p{r) we shall replace the values of the coordinates 
corresponding to Ss by those of S, justified by the same 
argument that p(r) is expected to vary over length scales 
which are much larger than d, we shall formally apply the 
second boundary condition directly on the surface S of 
the particle, rather than on Ss- Hence, the BCs take the 
following forms: 

p(|r| -^oo)^0 (6a) 



-D [h ■ Vp(r) 



IrGi; 



i^B o-T(r) 



(6b) 



= : i^B o- 



1 , r e catalytic cap , 
, otherwise . 



Examples of steady-state distributions p{r) of the so- 
lute (in units of po = i^b fi?i/D), obtained from a direct 
numerical integration of Eqs. ([5|) and ([5]), are shown in Fig. 
[21 One notes that as expected (sec Subsec. 13. ip there are 
significant density variations over the part of the surface 
covered by the catalyst. These are more pronounced (in 
addition to the larger amplitude of the distribution) for 
the "flatter" shapes, thus intuitively suggesting a larger 
resulting phoretic velocity for the more flatter shapes. In 
; the following subsections we derive the solute distribu- 
tions analytically for any spheroidal shape, which allows 
us to discuss quantitatively the effect of the shape on the 
phoretic velocity. 



3.3.1 Prolate particles 

For a prolate shaped {sr < 1) particle with the catalyst 
distributed over a cap-like region centered at one of its 
poles (see Fig. [1]), Eq. ^ subject to the BCs given by 
Eq. ^ is most conveniently solved in terms of the prolate 
spheroidal coordinates {£,,(f>,ri) (see Ref. [571 (a-)])- 



X = Ky^ (^2 _ _ ^2-) C0S(/) , 

y = K\/(^^ - 1)(1 - v'^)sm(j), 
z = kS^i] . 



(7) 
(8) 
(9) 



Z?V^p(r) = 0, r e outer region . 



(5) 



In these equations ^ > 1 and —1 < i] < 1 parameterize 
confocal ellipsoids and hyperboloids of revolution, respec- 
tively, which have their foci placed on the z axis a distance 
2k apart and symmetrically with respect to the origin O, 
while < < 27r is the azimuthal angle. The xz planar 
cuts through the iso-surfaces in terms of prolate spheroidal 
coordinates are shown in Fig. [IJb). The choice 



This equation is to be solved subject to the BCs of (i) zero 
density far away from the outer edge Ss of the surface 



Rl 



„2 



(10) 
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Fig. 3. Steady-state distributions p(r) of the solute (in units 
of po = vb <y Ri/D) around (a) a prolate {sr — 0.3), (b) a 
spherical (s,. — 1), and (c) an oblate (sr = 3.0) particle par- 
tially covered by catalyst {sh ~ 0.3). These distributions are 
obtained by numerical integration of Eqs. ([5]) and (|6]). In each 
case the color coding of the density ranges from blue (low den- 
sity) to red (high density) in accordance with the corresponding 
scale shown on the top left. Note the differences in the scales 
for (a), (b), and (c). 



ensures that the family of ellipsoids includes the one, 



which corresponds to the surface of the particle [^0 follows 
from Eqs. ©-((lOl) and z{^o, 77 = ±1; Kp) = The outer 
region thus corresponds to ^ > Co- Since the intersection 
of any hyperboloid rj = const from the family defined by 
Kp with the surface of the particle is a circle parallel to 
the equatorial xy plane, the cap-like region on E covered 
by the catalyst is parameterized by (^ = Co,0 < < 
27r, — 1 < ?7 < tjq), where 



- 1 < ?7o = -1 + s/i < 1 



(12) 

-Ri+ 



[which follows from Eqs. ([7|)- ([TT1) and z{(,q, rjo; Kp) 
h, see Fig. HJa)]. 

In terms of prolate spheroidal coordinates the solution 
p{^,r]) of Eq. which has no dependence on due to 
the azimuthal symmetry of the system, is finite at 77 = ±1 
(i.e., on the z axis), and satisfies the BC in Eq. ([Sa]). can 
be written as 



(13) 



£>0 



Pi and Qi 
Qeiw) 



qf^ with [SiEO] 

(£ + 2s)!(£ + s)! 



^ s!(2^ 



2sH 



2£- 1 



2^-5 



Pt-z{w) 



w e 



-2s-l 

Pe-i{w) 
\[- 



(14) 



1,11 



are the Legendre polynomial and the zeroth-order associ- 
ated Legendre function of the second kind of degree £ (Ref. 
[27l (b)]), respectively. The terms with Pf<o are, by defi- 
nition, identically zero. The coefficients ci are determined 
by the BC in Eq. (|6b|) . By noting that n = e^, where is 
the unit vector corresponding to the ^ direction, Eq. (j6b|) 
can be re-written as 



D 



= iyB(yT{ri;r]a) , (15) 



where /i^ is the scale factor corresponding to (with sim- 
ilarly defined scale factors /i^ and /i^) (Ref. (a)]): 




(16a) 



(16b) 



h^ = h^i^,r^) = Kpy/ie-m-V^), (16c) 

and the parametric dependence of the characteristic func- 
tion T(r) on 770 [Eq. ()6b|) ] is indicated explicitly. By using 
the orthogonality of the Legendre polynomials. 



Co 



> 1 



(11) 



2n + l 



(17) 
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ifherc Snm is the Kronecker delta symbol (Rcf. 23 (b)]), 



and by combining Eqs. (|T3| - (jl6ap . one finds for the coef- 
ficients cg: 



2 D Q',i^o)^oV^ 



(18) 



where Q^(0 = dQe{0/d^ and 



-1 

»?0 



dv^Ja-v'Piiv) (19) 
so that the coefficients q(^Oi'7o) arc dimensionlcss. 
3.3.2 Oblate particles 

For an oblate shaped {sr > 1) particle with the catalyst 
distributed over a cup-like region centered at one of its 
poles [see Figs. [1] and [5], Eq. ([S]) subject to the BCs given 
by Eq. ([6]) is most conveniently solved in terms of the 
oblate spheroidal coordinates (C) x) (see Ref. [171 (a)]): 



^ = «;V(C' + I)(l-X')cos0, 
y = nViC + 1)(1 -X^)sin(/), 
z = kCx ■ 



(20) 
(21) 
(22) 



In these equations C ^ and — 1 < X ^ 1 parame- 
terize confocal ellipsoids and half-hyperboloids of revolu- 
tion, respectively, which have foci placed symmetrically 
with respect to the origin O on the x axis a distance 2 k 
apart, while < < 27r is the azimuthal angle. Cuts by 
the xz plane through the iso-surfaces in terms of oblate 
spheroidal coordinates are shown in Fig. [Ijc). The choice 



Rl- Rl = Ri^Jsl-l 



(23) 



ensures that the family of ellipsoids includes the one, 

1 



Co 



>o, 



(24) 



which corresponds to the surface of the particle; Co follows 
from Eqs. (HOll-dMl) and z(Co, X ±1; ^o) = Ri- The outer 
region thus corresponds to C > Co • Since the intersection of 
any hyperboloid x — const from the family characterized 
by Ko with the surface of the particle is a circle parallel to 
the equatorial x y plane, the cap-like region on S covered 
by the catalyst is parameterized by (C = CoiO < < 
27r, -1 < X < Xo), where 



Xo 



-l + Sh- 



(25) 



This follows from Eqs. (gOll-dMl) and z(Co, Xo; ^o) = -Ri + 
h [see Fig. ma)]. 

In oblate spheroidal coordinates the solution p(C, x) 
of Eq. ([5]), which has no dependence on due to the az- 
imuthal symmetry of the system, is finite at x = ±1 (i.e., 
on the z axis), and satisfies the BC in Eq. (|5a)) . can be 
written as [29l 



(26) 



l>0 



where i = and [see Eq. ([H))] 

(^ + 2s)!(^ + s)! 
s!(2i' + 2.s + 1)! 



Q,(iC) := 



s>0 

i (arccot C)P£(iC) 
'3(731)^^-3(10 ■ 



2£- 1 



^^-i(iC) 



(27) 



The second equation, for which we have used (1/2) ln[(2 + 
l)/(z— 1)] = arccoth(z) and arccoth(iz) = — iarccot(z), 
emphasizes that QiiiQ is well defined for all C S K |29j . 
This is needed because, for Sr > V2, Co is smaller than 1 
[Eq. ([M)) ] and thus the point z = i, where the Legendre 
differential equation is singular, lies inside the domain of 
the solution (in contrast to the case of a prolate, for which 
5o > !)• The coefficients are determined by the BC in 
Eq. ()6b|) . Note that because the density p{(, x) G while 
according to Eq. ([?7)) Qe{K) is real (imaginary) for £ odd 
(even), the coefficients nig in the series representation, Eq. 
((26)) are real (imaginary) for £ odd (even). 

By using n = e^, where e^; is the unit vector corre- 



sponding to the C direction, Eq. 
as 



6bl) can be re-written 



D 



hdCo,x) V dC 



where /ij , h^, and /i^ with 



= j.b(tT(x;xo), (28) 



C=Co 



= hiC^x) = 




= hdc,x) = KoVic + m-x'), 



(29a) 

(29b) 
(29c) 



are the scale factors corresponding to the e^, e^, and 
directions, respectively (Ref. [571 (a)])- By using the or- 
thogonality of the Legendre polynomials [Eq. (fT7)) ]. and 
by combining Eqs. ([55 |) - (|29a[) . one finds the following ex- 
pression for the coefficients mg: 



^KCcXo) 



vb^tRi i2£+l)eg{Co,Xa) 
D 2CoVC^[9cQ^(iC)]c=Co 



— — "^HCo,Xo) , 



(30) 
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where 



3.5 Phoretic velocity of an oblate object. 



-1 



dx ^Q+X^Pe{x) (31) ment is dS = dx 

so that the coefEcients TO£(CoiXo) are dimensionless. 



Similar to the calculation in Subsec. 13.41 for an oblate 

shape one has [3T] : 

\ d 1 ^ 

(i) = - — — h eA - — — — and the surface area clc- 

^ h^dx ocj) 



3.4 Phoretic velocity of a prolate object. 



, , xdrx + ydry + zdrz 
(n) n • r = • r = — ^ H ^ = — , and 



C 



(iii) • V^p(C, X) = {e. ■ 6x) = ^ ^xP(C, X) . 



By using the expression in Eq. (|4]) for the phoretic velocity 
of an o61ate object one obtains: 



From differential geometry one has |3T| : 

Id Id 
(i) = e„ — — — h — — — and the surface area ele- 
hri orj 84) 

ment is dE = h^i drj dcp because the plane tangent to 

the surface of the particle is spanned by the unit vectors 

and of the t] and </> directions, 



(ii) n • r = • r 



xdfx + ydev + zdez kIS, 
— ^ ^ ^ — -r^, and 



(iii) e. • V^p(e, ^) = (e, . e,) Mil^ = M Q.pi^, r^) . 

By using the expression in Eq. Q for the phoretic velocity 
of a prolate object one obtains: 



Vr 



Ri J (Co - t) 

-1 

-21/0^2(^0-1) E [^2^+1(^0,^70) 



t>0 



{21+1 



(Co) / drj 



■P2i+i{v)]- (32) 



i-x' 



-1 



Ho'(Co^ + l)' 



dx ■ 



-1 



p(Co,x) 



-21/oCo'(Co+l) E[™2£+i(Co,Xo) 



X Q2£+i(iCo) / dxj^Tf^P2e+i{x)]- (34) 



After replacing rhi, Co, and xo by the corresponding ex- 
pressions in Eqs. (pO)) . (|24p . and ((25|) . one obtains as final 
result the phoretic velocity Vob = Vob{sr, Sh',Vo). Similar 
to the case of a prolate object, the terms with even i do 
not contribute to the phoretic velocity, i.e., mi even = 
because Pi{r]) is an even (odd) function of x for ^ even 
(odd). Note that the expression in Eq. can be ob- 
tained from Eq. p2[) by the mapping ^ i^; this is in 
agreement with similar observations regarding solutions of 
the Laplace equation in spheroidal coordinates (see, e.g., 
Ref. [32]), which is a welcome consistency check for our 
results. 



where 



Fo = 



D 



(33) 



defines the velocity scale. Vb is expected to be of the order 
of fim/s, but because of its dependence on b it is difficult 
to provide a theoretical estimate for it. The second equal- 
ity follows from an integration by parts, while the third 
one uses the series expansion in Eq. for p{£,o,r]), Eqs. 
pH]) and pT|) to replace Kp and s^, respectively, and the 
fact that Pe{ri) is an even (odd) function of 77 for £ even 
(odd). After replacing ci, ^q, and 770 by the corresponding 
expression in Eqs. ([T5|) . pT|) . and P^ . one obtains as fi- 
nal result the phoretic velocity Vpr = T^r(sr, s^,; Vb) as a 
function of the geometrical parameters Sr and Sh, as well 
as of the velocity scale Vq. 



4 Discussion 

We first note that the length scale Ri does not enter ex- 
plicitly into the final expression for the velocity [Eqs. 
and (j34p ] which shows that within the assumptions of the 
model the phoretic velocity of objects with the same as- 
pect ratio but different linear sizes is the same (in agree- 
ment with the conclusions of Ref. [11] for spheres and 
cylindrical rods) . 

The second observation concerns a symmetry with re- 
spect to the area covered by the catalyst. The velocity 
depends on the ratio Sh, or equivalently rjo or xo, via the 
integral in Eq. (fT9|) and Eq. ([3T|) . respectively, which enters 
into the coefficients of the corresponding series expansion 
of the density. With Pei-rj) = -Pe{v) for £ odd Eq. ^ 
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yields for < w < 1 and I odd 



-1 

w 
-1 

= = w) , (35) 

1 

because / d-q \/Q--rfPe{i]) = for ^ odd. A similar re- 

-1 

lation holds for the case of an oblate shape: 

fflCcXo = -u') = q(Co,Xo = u;) • (36) 

Since Sh = l + f]o = 1 + XOj this means that the coefficients 
in the series expansion in Eqs. ([5^ and and thus the 
velocity of the particle, are the same if the catalytic cov- 
erage is less or more than half of the particle by the same 
amount w. Equivalently, it means that if the inert and the 
catalytic characteristics of the two parts of the surface are 
interchanged, the velocity is the same but the direction 
of motion is reversed. Therefore, as a function of s/j, Vpr 
and Vob have an extremum at s?i = 1 (?7o = Xo = 0), i.e., 
if half of the particle surface is covered by the catalyst. 
For Sfi = (no catalyst) or s;i = 2 (entire particle surface 
covered by the catalyst) and a finite number density a of 
catalytic sites, the coefficients C21+1 and m2i+i are iden- 
tically zero and the velocities predicted by Eqs. ([5^ and 
vanish as expected because in these two cases there 
are no gradients of the density of the product molecules 
along the surface of the particle. 

The dependences on ^0 and r/o as well as on C,q and 
Xo (i-e., on s,. and s^) of the series representations of the 
velocity [Eqs. ([5^ and ((M)) ] are very complicated and in 
the general case we have not been able to obtain a solution 
in closed form. Therefore, we shall study numerically the 
cases of generic prolate (0 < Sr < 1) and oblate (1 < 

< 00) shapes, and we shall complement the analysis 
with analytical results for the limiting cases of a spherical 
{sr = 1), a needle-like (s^ 0), and a disk-like (s^ — > 00) 
shaped particle. 

For given ^0 and r/o (Co and xo) we approximate the 
velocity Vpr {Vob) by keeping terms up to ^ = 15 in Eq. 
((32|) [Eq. (p4)) ]. This provides a good approximation for 
aU values (a) 1.001 < ^0 < 80 (i.e., 0.0451 < < 0.9999) 
and — 1 < ?7o < 1 (prolate shapes) and (b) 0.01 < Co < 80 
(i.e., 1.00008 < Sr < 100) and -1 < xo < 1 (oblate 
shapes) which we have tested (in the sense that including 
five additional terms leads to changes in the value of the 
velocity smaller than 10~^) [33]. The results for V^r/Vb 
and Vob/Vb are shown in Fig. S] as functions of ryo and xo 
for several values of ^0 and Coj respectively. The data span 
the whole range of interest, from ^ 1 (i.e., ^0 ^ 1), 
corresponding to a "cigar" shaped rod, to ^ 1 (i.e., 
Co ^ 1), corresponding to a disk, through ~ 1 (i.e., 
^0 ^ 1 or Co ^ 1) corresponding to a slightly deformed 
sphere. 




-1 -0.5 0.5 1 

Fig. 4. Dependence of the scaled phoretic velocity [Eqs. p2|) . 
(|34p . and (|33|l ] on the fraction = 1 -I- t;o = 1 -I- Xo of the 
surface of the particle covered by the catalyst (t^o, Xo = — 1, 0, 1 
correspond to a particle surface with no catalyst, with the lower 
half covered by the catalyst, and completely covered by the 
catalyst, respectively), for (a) a prolate-shaped particle with 
aspect ratio parameter ^0 = (1 - s^)"^''^ = 80 (sr = 0.9999), 
2 {sr = 0.866), 1.1 {sr = 0.42), 1.01 (s^ = 0.14), and 1.005 
(sr = 0.099); and (b) an oblate-shaped particle with aspect 
ratio parameter Co = (-1 + Sr)"^''^ = 80 {sr = 1.00008), 1.0 
{sr = 1.41), 0.5 {sr = 2.24), 0.25 {sr = 4.12), and 0.06 [sr = 
16.7). The numerical results are shown as symbols; the dashed 
lines connecting the symbols are guides for the eye. Also shown 
(solid line) is the exact result [Eq. (|65p ] corresponding to a 
spherical particle (^0, Co — > 00). 



Figure m shows that the dependence of V/Vq on 770 (or 
Xo) is qualitatively the same for all values of Co and Co- 
The ratio V/Vq is always negative, which means that V 
and Vo have opposite signs. Therefore, for repulsive ef- 
fective interactions between the particle and the product 
molecules, i.e., in the case that 6 < and thus Vb < 0, 1/ 
is positive and the particle moves in the positive z direc- 
tion, i.e., "away" from the catalyst covered end. Similarly, 
for attractive effective interactions the motion will be di- 
rected "towards" the catalyst covered end. As expected 
the velocity curves are symmetric with respect to 770 = 
and Xo — (i-C-, s/i — 1), which is a highly welcome 
check for the numerical calculations. For a prolate, the 
minimum at 770 = slowly rises towards zero upon de- 
creasing Co towards Co = , so that the velocity vanishes 

for Sr = 1 — C(7^ 0; as discussed below, this is in 
qualitative agreement with Refs. pilSllllj. In contrast, for 
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an oblate the absolute value of the velocity increases upon 
decreasing Co (at fixed xo) towards a non-zero and finite 
limiting value which corresponds to that for an infinitely 
thin disk (oblate of equatorial diameter much larger than 
the polar one) [see, c.f., Subsec. H^ . (Note that the veloc- 
ity of such an object remains finite because the product 
molecules can still "escape" to both sides of the object 
due to the assumption of an unbounded solvent.) There- 
fore, we conclude that for a given coverage by catalyst a 
flatter shape of the particle leads to an increased phoretic 
velocity. This can be understood intuitively by analyzing 
the case of a single, point-like catalyst source placed at 
one of the poles of the particle. The spherically symmet- 
ric product number density distribution due to a point-like 
source but in the absence of the particle will be deformed 
over a larger extent by an impenetrable boundary which 
is more elongated transversal to the radial direction than 
by one more elongated along the radial direction, and thus 
the gradient along the surface of an oblate particle is ex- 
pected to be larger than the one corresponding to a prolate 
particle. Finally, we note that (i) for large values of or 
Co the numerical solutions are in perfect agreement with 
the known analytical result [Eq. (j65p ] for a sphere [TT] : 
and that (ii) for very elongated objects (rods or oblates) 
the results should be considered only as a first-order ap- 
proximation. This is the case because the analysis leading 
to the expression in Eq. ([T|) approximates the fiow within 
the thin surface film as locally planar; such an assump- 
tion may break down near the poles for prolates or near 
the equator for oblates. There for very elongated particles 
(sr- <C 1 or Sr ^ 1) the curvature of the surface of the 
particle is very high and cannot be neglected a priori. 

Further insight can be obtained from an analytical 
study of the behavior in the limiting cases of a slightly 
deformed sphere: ^o ^ 1 {sr ^ 1) or Co 3> 1 (s,. > 1), of a 
needle-like shape: Co 1 {sr ^ 1), and of a disk: Co ^ 1 
{sr 3> 1) which we shall discuss now. 



4.1 Limit of a slightly deformed sphere 



Accordingly, Eq. reduces to 



Vpri^O > 1) 
Vn 



-2/1 ^-2n1/2 71 (Co)Ql(^o) 



Mo) 



Furthermore, Eq. (jl9|) implies 



(37) 



71 (Co » 1) =Co 



+ § Co"' j dv rf r{rr, Vo)Pi{v) + 0(^0^") 
-1 



(38) 



note that the first integral on the right hand side is equal 
to twice the velocity (in units of the characteristic ve- 
locity Vb) corresponding to a spherical particle [Eq. ((65|) . 
Appendix IB|) ]. Due to 



l(Co) 



CoQUCo) 



= -2 + Y^Co"' + ^(Co"') (39) 



and Vsph/Vo = - l)/4 one finds with r]^ < 1 [Eq. (fH)) ] 



Vpriio » 1) 
Vo 

Vs 



5if - 6 



Vo 20 

-1 

Vsph (-^ _ 17 - bril ^ Vsph 



20 



^o]+0{i^^) 



20 



Vo ' 



(40) 



We first consider the case in which the limit towards a 
sphere is taken through prolate shapes, i.e., Sr — > 1^. 
Noting that in the limit of a sphere Sr = R2/R1 — >■ 1~, 
which is equivalent to Co — ^ co [see Eq. (TTlT) ]. the integral 
in the last line of Eq. ([5^ can be approximated as: 



(C(f - v^y 



.-4 

so 



^0^ 



dririP2i+i{Ti) + 0{i^^) 



-<5,.o + 0(Co-^) 



where the last inequality follows because the velocity ra- 
tios are negative quantities. Equation (j40p thus implies 
that the absolute value of the velocity of a sphere slightly 
deformed towards a prolate shape (Co >• 1) is smaller than 
that of a sphere and becomes equal to it in the limit of a 
vanishing deformation (Co — >■ 00), in agreement with the 
numerical findings shown in Fig. Sl^a). 

A similar calculation can be performed for the case 
in which the limit is taken through oblate shapes, i.e., 
Sr — > 1^, which is equivalent to Co — >■ 00 [see Eq. ((24| ]. 
by starting from the integral in the last line of Eq. 
Following the same steps as above and noting that 



Qi(iCo) 



Co[5cQi(iC)]co 



Co»l 



= -2 + I^Co"' + ^(Co"') (41) 
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we obtain with Xo — 1 due to xo = Vo [see Eqs. (fT^ and 



(Hi 

VobiCo » 1) 
Vo 



20 



Vo 



+ 0(Co"') 



1 - 



20 



Vn 



20 



^0 ' 



(42) 



because the velocity ratios arc negative quantities. Thus 
the absolute value of the velocity of a sphere slightly de- 
formed towards an oblate shape (Co 3> 1) is larger than 
that of a sphere, while in the limit of a vanishing defor- 
mation (^0 — >■ oo) it indeed reduces to the velocity of a 
spherical particle, again in agreement with the numerical 
findings shown in Fig. lUJb). 



4.2 Limit of a needle-like particle 

In this subsection we focus on estimating the asymptotic 
behavior of the velocity of a half-covered (sh = 1, ?7o = 
0) prolate particle with increasing elongation towards a 
needle-like shape (s^ <C 1, i.e., ^ !)■ As shown in Fig. 
IDJa), for given < 1 (or, equivalently, > 1) the abso- 
lute value of the velocity for any value ijo 7^ is smaller 
than the one for 770 = 0. Thus the asymptotic behav- 
ior of |lpr(^o %) I has an upper bound given by 
lV(eo^l+,r/o = 0)|. 

Eq. ([5^ may be re- written as: 



1 

VpriCo,m = o) = ^Jdv 



77'= 



= ^ [PUo,?7 = 1) - /3(?o,?7 = -1)] 



^2 ) 9nP{Co,v) 



2i?i 



ieo 1) 



£ odd 



X! C£(?o,??o = 0)Q£(Co) dr]-p^ 

5] Q(eo,??o = 0)Q,(^o) 
1 dP, 



1 m 

— dr] 



e odd 



1 - ^ / dr, 



^0 - V'^ dri 



(43) 



where we have used Pi{l) = 1 [23 (c)] and the fact that 
for even indices (. the integrals are identically zero because 
in this case PeXv) is even function of r] so that the 
integrands are odd functions of r/. For odd i the derivative 



dP( / drj is a polynomial of order £ — 1 containing only even 
powers of 77. Thus it can be written as 



Pi-iiv) := 



dPj 
drj 



mo) + ieo-v')qi-3im^o), (44) 



where 



mo) 



dPe 
drj 



(45) 



By construction the second term in Eq. (|44)) must vanish 
for 7? = ^0- Since dPi/dr/ depends on r] only via 77^, this 
vanishing must exhibit a prefactor — ry^ multiplying a 
polynomial qe-3{r];^o) of degree (£ - 3) with qm<o = 0. 
Accordingly, the integral in the last equation of Eq. 
is approximated by 



drj-j- 



— = ? — In c , -I + / a'7 9f-3 '7;Co) 

'7 a?? Co Co + 1 7 

-1 

-/3,(eo = l)ln(Co-l). (46) 



This implies that in the limit Co ~^ l"*" the second term in 
the square bracket in Eq. (j43p vanishes ~ (Co — 1) ln(Co — 
1) and therefore it is a subdominant contribution to the 
quantity in the square bracket. Thus Eqs. and (HH) 
lead to the following approximation for the velocity of a 
very elongated prolate towards a needle-like shape (Co 
1+): 

^pr(C0> 1,^70 = 0)^^ 

X ^^(^"'^" = O)Qf(Co) [1 + A(l) (Co - l)ln(Co - 1) 

i odd 



^ E c.(Co,'7o = 0)Q,(Co) 



i odd 
Eq. Ul Vb 



I odd 



(2^+1) 



<3«(C0) „^ 

7777T7KCo,f7o = 0) 



(47) 



Within this approximation the velocity of a very thin, 
needle-like prolate object is proportional to the difference 
between the values of the product molecules density at 
the two ends (77 = ±1) of the object. (Note that because 
of this subtraction only those terms with odd I occur in 
the series representation of the velocity.) This is similar 
to the results for a cylindrical thin rod postulated in Ref. 
[2] and derived in Ref. [11] by invoking a "slender-body" 
approximation. 

Although the above approximation leads to a signifi- 
cantly simplified expression for the velocity, we have been 
unable to further simplify the resulting series and thus we 
have studied it numerically. The conclusion of this analy- 
sis, the details of which are presented in Appendix [Cj is 
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that for Co Jc; 1 the series in Eq. (|T7)) behaves as 



e odd 

^ (Co-l)ln(eo-l)/(eo) 



mo) 



where 



/(Co>l)- -1.5X [-ln(eo-l)]- 



(48) 



(49) 



Eqsuations (|47)) and (|49)) therefore render in the limit 

x[-ln(Co-l)]-"-'<0. (50) 

In agreement with the behavior observed for a general pro- 
late shape [see Fig. Ufa)], the velocity (in units of Vb) is 
negative. For ~^ 1^: the velocity vanishes faster than 
VCo — 1 lii(Co — l)j i-e., more rapidly than the behavior 
predicted by Refs. [TTll^. This difference can be either due 
to an intrinsic difference between the cylinder-like and the 
needle-like shapes, or, most likely, due to the additional 
approximations employed in Refs. [TTJ[2] upon computing 
the solute density distribution (such as using a distribu- 
tion of point sources in unbounded space and the absence 
of sources on the flat ends of the cylinder) . 

We note that in terms of an effective power law the 
value 0.9 of the exponent in Eq. ([^^ provides a very good 
approximation for the behavior of /(^o) over the physi- 
cally accessible range of values ~^ l^j which can be es- 
timated to be bounded from below by — 1 + 10~^ [i-e., 
the value corresponding to a 10 /xm long carbon nanotube 
of 1 nm diameter {sr = 10~^)]. However, in a strictly 
mathematical sense, the limiting behavior of /(^o l"*") 
appears to be not given by Eq. (|49|) because one finds that 
the exponent decreases as the range under considera- 
tion corresponds to smaller and smaller values; e.g., the 
expoenent reaches the value 0.7 for <^o Jc] 1 + 10^^^*. 



4.3 Limit of a disk 

In this subsection we focus on the case of a half-covered 
{sh = 1, Xo = 0), very flat oblate (disk-like shape, Co *C 1) 
because among the class of spheroids wc have studied this 
is the one which exhibits the largest absolute value of the 
velocity. 

For Co ^ 1; the coefficients e„ with an odd index n = 
2£ + 1 [Eq. ([5T|) ] can be approximated by 



xo=o 



e2£+i(Co < l,Xo = 0) 



= -\jdxPi{x)P2i+i{x) = -\5ifi. (51) 



(According to Subscc. [331 the coefficients e2i do not con- 
tribute to Vob-) Therefore in this limiting case only the 
first term in the series representation given in Eq. 
contributes. Since 



Oi(iCo) 



[9cQi(iC)]c=Co 



-1 -I- Coarccot(Co) 



Co 



1 + 
2 



arccot(Co) 



Co«l 



(52) 



and 



dx 



■Piix) 



- Co ^arccot(Co) , (53) 



one obtains 



K&(Co < l,Xo = 1) - - 



Vo 



— !• 



Vo for Co ^ . 



(54) 



Thus in our context the maximal velocity of a spheroidal 
object is |Vb/2| and it is realized for a particle with a disk- 
like shape and such that one side is covered by catalyst 
and the other one is inert. This is in agreement with the 
numerical results in Fig. |31 

The general case Xo 7^ can be studied in a similar 
way, but due to the fact that in this case all coefficients 
££ with £ odd will contribute [see Eq. ([5T|) with xo 7^ 0] to 
the series representation in Eq. (j34[) . the resulting velocity 
has to be calculated numerically. 



5 Summary 

We have studied the diffusio-phoretic velocity of a spheroidal- 
shaped particle (Fig.[T|) which self-propels by creating gra- 
dients of product molecules in a surrounding, unbounded 
solvent (Fig. [3]). Our calculations yield the dependence of 
the velocity on the shape of the particle, i.e., the aspect ra- 
tio between the polar and the equatorial diameters, and on 
the fraction of the surface of the dissolved particle which 
is catalytically active in providing product molecules via 
chemical reactions. 

The analysis of our model is based on recasting it in the 
framework of the standard theory of phoresis. In this con- 
text we have critically analyzed the assumptions involved 
in such an approach, some of them already present in the 
standard theory, others arising as a result of the mapping 
of such "active" surface particles into the framework of a 
theory developed to describe the case of inert particles im- 
mersed in a pre-defined, externally controlled concentra- 
tion gradient. Within the confines of the standard theory 
of phoresis, we have shown that, irrespective of the shape 
and of the fraction of the surface covered by the catalyst, 
the phoretic velocity of the particle depends on its geom- 
etry only via the aspect ratio but it is independent of its 
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absolute size. The numerical analysis of the series repre- 
sentation for the phoretie velocity has been complemented 
by analytical results for the asymptotic cases of spherical 
and needle-like particles (Fig. [3]). For a given coverage of 
the particle surface by the catalyst, the absolute value of 
the velocity of a prolate particle is maximal for an almost 
spherical particle and decreases gradually towards zero 
with increasing elongation of the prolate towards a needle- 
like shape [Fig. IH^a)]. We performed thorough numerical 
studies (Fig. [5]) of the decay of the velocity with increas- 
ing elongation {Ri ^ R2, see Fig. [T]) and concluded that 
over the whole physically accessible range Co ^ 1 + 10^*^ 



-In^ 

Ri 



it decays effectively ~ [ — In — 
\Ri Ri^ 

R f R \ ^'^ 
— — - — In — - ) , which is faster than the previously 

^1 V Ri J 

f R2 R2\ 

predicted — In — behavior pTll^lB] ; this difference is 
\Ri RiJ 

most likely due to the approximations employed in these 
previous calculations. Our numerical study also indicates 
that this conclusion of a faster decay than previously pre- 
dicted holds in the limit (of pure mathematical interest) 
,^0 — >■ l"*" (see the inset in Fig. [S]), but we could not de- 
termine the exact analytical form of this decay. In con- 
trast, an oblate-shaped particle moves faster with increas- 
ing flattening from an almost spherical towards a disk-like 
shape [Fig.lU^b)]. For a given shape the maximal absolute 
value of the velocity V in units of the characteristic veloc- 
ity Vq is always attained at half-coverage by the catalyst 
and varies from (needle) over 1/4 (sphere) to 1/2 (disk). 
Therefore, experimental realizations of such self-propelled 
particles call for a compromise on one hand between the 
increased speed of flatter particles and their decreased uni- 
directionality due to, e.g., the unavoidable thermal noise 
of the solvent and the density fluctuations of the product 
molecules, and on the other hand between the increased 
stability against rotations of more elongated particles and 
their reduced velocity. The results for the phoretic velocity 
V/Vq shown in Fig. 2] allow one to directly compare them 
with experimental realizations, e.g., confocal microscopy 
studies of micron sized needle-, sphere-, and disk-shaped 
particles, both qualitatively - the symmetry with respect 
to the half-covered case and the dependence on the aspect 
ratio Sr — R2/R1 of the shape - as well as quantitatively: 
the velocity of a half-covered disk is twice that of a spher- 
ical particle of the same radius. 



ity and the uni-directionality of the motion for spheroidal 
particles with axial symmetry. 
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A Derivation of the phoretic velocity 

The flow field u(r) in the outer region is the solution of 
the incompressible, force free Stokes equations 



v-n = o, V-u = 0, 

subject to the boundary conditions 



= V + v, 



(55) 
(56) 



n :~ —pI+iiS is the corresponding pressure tensor, where 
p is the hydrostatic pressure and S is the shear stress 
tensor, i.e., Sa/s = dua/dxp + du^jdxa- Owing to the 
linearity of the Stokes equations, we can write the solution 
as u = u' + u" and n = n' + n", where u' and u" are 
the solutions of the incompressible Stokes equations which 
vanish at infinity and satisfy the boundary conditions 



u 



V, and u" 



= V. 



(57) 



respectively, while 11' and fl" are the corresponding pres- 
sure tensors. By using the condition that the motion of 
the particle and the outer hydrodynamic fiow are such 
that there is no net force acting on the system composed 
of the particle plus the surface film, and by replacing n" 
with the expression for the pressure tensor on the surface 
of an ellipsoid translating through an unbounded fiuid at 
rest [Ml, 



3 V„ ^ 



(58) 



As we pointed out in our discussion of the connection 
between the phoretic slip and the number density gradi- 
ents of product molecules (Subsec. 13. 2[) . a natural exten- 
sion of the present study would be to consider in detail 
the generic case in which the effective interaction poten- 
tial between the product molecules and the particle differs 
on the catalyst-covered part of the particle surface from 
that on the bare one. Other questions of further interest 
are the influence of curvature on the phoretic velocity in 
the case of very elongated or very flat spheroidal particles 
and the influence of external boundaries on both the veloc- 



one obtains 



di:fin' 



ds n n" = k(*) V 



(59) 



n is the unit vector of the direction normal to Ss (oriented 
towards the fluid) and K*^*' is a constant diagonal tensor 
which depends on the diameters of the particle only; its 
explicit expression is not needed in the following. 

On the other hand, because the pairs (u' , 11') and 
(u" ,n") are, by construction, solutions of the force-free. 
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incompressible Stokes equations and decay at infinity, they 
satisfy Brenner's (or Lorentz's) reciprocal theorem [181 
(b)]: 



respectively, one finds: 



dSnW u" 



dSnll" u'. 



(60) 



By using Eq. dSS]) to replace fl" and the BCs in Eq. ([57]) 
to replace u' and u", and by noting that V and K*-*-* are 
constant with respect to the integrations, one obtains 



V 



dSnll' ^---^K'-'^Y ffdS(h-r)w, 



3 V 



(61) 



By using Eq. ((59)) for the rhs of Eq. (1611) and then replacing 
in the calculations Us by one obtains Eq. (U) for the 
phoretic velocity. 



B Phoretic velocity of a spherical Janus 
particle 

The case of a spherical particle (s^ = 1) with a catalytic 
cap centered at a point chosen as one of its poles (see 
Fig. [U has been discussed in Rcf . [TTj . Thus we provide 
here only a brief summary of the results for reasons of 
completeness, further referencing, and provision of some 
remarks. Equation ([S|) subject to the BCs given in Eq. 
(O can be conveniently solved using the standard polar 
spherical coordinates (r, 0, </>). The solution p{r,9) of Eq. 
(O (there is no dependence on cj) because the system has 
azimuthal symmetry) which satisfies the BC in Eq. (j6a[) 
is 

where the coefficients di are determined by the BC in Eq. 
(j6bp . Since the cap- like region on U covered by the catalyst 
can be parameterized hy {r = Ri,9q < 9 < n,0 < (f) < 
2tt), where 

cos 00 = -1 + s/i , (63) 

one obtains 



7T 

Vsph ^ j d9 sin 9 cos 9 p{Ri , 



2,.T, ^ (l--s,0^-l ^^ 
-Vodi{sh) = Vo . 



(65) 



In Eq. ([55)) the first equality has been derived in Refs. [TJ 
111) . The second equality follows by using the series expan- 
sion in Eq. ([52) for p{Ri,9) and due to cos 6* = Pi (cos 6*) 
so that in the expansion all terms with i ^ 1 vanish. We 
note that, as for the general prolate or oblate spheroid, 
there is no explicit dependence on Ri [TT]. Equation ([55)) 
provides the following conclusions: 

(i) Since < s/i < 2, Vsph and Vq have opposite signs. 
Thus for repulsive interactions between the particle and 
the product molecules the motion will be in the positive 
z-direction because b and thus Vo are negative [T^[7) ; sim- 
ilarly, for attractive interactions, the motion will be in the 
negative z-direction (see Fig.[T]). 

(ii) As expected on basis of symmetry considerations, the 
maximum velocity VJ^^^ = — Vo/4 occurs for Sh = 1, i.e., 
if the catalyst covers just a hemisphere. 

(iii) For any finite value a of the density of catalytic sites 
the phoretic velocity vanishes if the surface area covered 
by the catalyst tends to zero, i.e., Vsph{sh — > 0; cr < oo) — > 
0, as expected intuitively. 

(iv) Since 2'KR\sh<J is the total number of catalytic sites 
in the area covered by the catalyst, the case of a sphere 
with a single catalytic source corresponds to the limit 
{sh — > 0, (T — > oo} such that 27r i?fs/i(T = 1. The phoretic 
velocity in this case is 



Vsph{{sh -> 0, CT 0O,2TTR\shU = 1}) = - 



in agreement with Ref. [7]. 

(v) If a general axially symmetric distribution of catalytic 
activities i^Bf{9) is considered, as in Rcf. [TT], the 

second equality in Eq. ([65)) still holds with the coefficients 
di redefined by the particular choice of the distribution 
considered by replacing T{9) i— )■ f{9) T(9) in the first equa- 
tion in Eq. ([M)). 



2i+l R{ 



+2 



2(£-|-l) D 



w 

D 



i+2 



d9 (sin 61) UB cr T{9, 9q) Pe (cos 9) 



2e + i 

2{l+l) 



■ de{sh 



+1 



dxT{x, cos 9o) Pe (x) 



C Numerical analysis of the velocity in the 
limiting case of a needle-like particle 



According to Eq. ([47| . the behavior of the velocity of a 
(64) prolate object in the limit ^ 1^ is determined by the 
series 



Withn-r = rande^-V^p(i?i,6') = -{siii9/ Ri) d0p{Ri,9), 

and by using the expressions in Eqs. and ([6b| for the f{£,o) {21 + 1) 

phoretic velocity and for the characteristic function T(r), 



l£i^o,Vo = 0) 



(66) 
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Because Qe{£,o 1^) ^ ln(Co — 1) [see Eq. (|Ti)) ]. one may 
expect that for ■^o ^ 1 the series above varies as 



Legcndrc polynomials with odd index £ in terms of the 
hypergeometric fmietion 2F1 [30] : 



/(eo^l+) ^ (eo-l)ln(Co-l) 



I odd 

=:(eo-l)ln(Co-l)/(eo-l), (67) 

where the series defining the prefactor /(I) = /(^o — 1) 
has to be calculated numerically. However, within the lim- 
its of numerical accuracy it turns out that /(I) = so 
that Eq. ([57)) does not capture the leading behavior of 
/(^o ^ 1^)- Most likely, the reason for the failure of 
this approximation is that the series representation of the 
product density [Eq. (|13p ] is not uniformly convergent 
(and actually completely breaks down) at = 1, where 
the differential equation of the Legendre functions is sin- 
gular and Qi{S,()) diverges; consequently, there is no war- 
ranty that the limit ^0 ~^ l"*" can be taken term by term. 
The result /(I) = thus suggests that instead the series 
has to be first summed up for general ^0 £md only then 
the sum can be evaluated in the limit ^0 — ^ l"*"- (See also 
Ref. [5], in which the thin rod limit of a vanishing ratio 
between the radius and the length of a cylinder [i.e., Eq. 
(4) therein] could be taken only after calculating the ve- 
locity as a function of this ratio [Eq. (3) therein]; similar 
arguments apply to the derivation of Eq. (16) from (14) 
in Ref. [H] ) 

We therefore proceed with a numerical study of the 
series /(^o) [Eq. ([66| ] as a function of ^0 ^ 1- By using the 
recursion relation satisfied by the Legendre functions 

m 



1 



[Q(+i{w) - wQi{w)\ 



(68) 



and the series representation [Eq. (jl4p ] for Qi{w), one 
concludes that for ^0 > 1 and f » 1 the behavior of the 
ratio involving the Legendre functions is given by 



mo) 



2^0 (e, 



1) 1 



£»1 



(1-2^2) i+1 



< 0. 



(69) 



Turning now to the behavior of Ii{^o) 
for Co ^ 1: we first note that [Eq. ([H 



mo)\ < / drj 







<^oJ dij\Pe{v)\<Co- 

-1 



(70) 



For large i, i.e., i > 50 and values of ^0 very close to 1, i.e., 
^Q — 1 < 10~^, an accurate direct numerical determination 
of Ie{£,o) [Eq. (|19p ] is extremely difficult mainly because 
of the oscillatory behavior of the Legendre polynomials, 
and one thus has to find a way to reformulate the inte- 
gral. To this end we make use of the representation of the 



Piirj) ^ (-1)^ 



X77 2F1 



2^-1 f LJ.\ 



at Peiv) , (71) 



where at denotes the prefactor and 

Pi{v) = ^7 2-^1 ( 2 ^' 2' '^^ j ^ polynomial of 

order £. The hypergeometric function 2-^1 (w) obeys the 
relation [27l (d)] 

d , 2\ 
— 2F1 [a,b;c;w j 2 — w 
dw c 

X2 Fi (a+ 1,6+ l;c+ l;w2) , (72) 

with 2i^i (0, 6; c ^ 0; w^) = 1 and 2^'i(a, b; c; 0) = 1, and it 
is well defined, as well as its derivatives, at w = 1 (being 
polynomials). Accordingly, for £ odd, Ii{^o) can be com- 
puted via successive integrations by parts: 



dri 



X 2F1 



1 £ 



-\^l + a,{£^ml~l?" 



'2 
+ I3 





X J drj 

-1 



£ 

^+2 



1 £ 
2^2 



£ ^ 1 £ 3 
X 2F1 { ^ + 1,- + 1 + 1;- + 1;77^ 



+ ai(£,0)(eo'- 1)3/2 + ai(£,l)(eo'- 1)5/2 
3/ UJ V2 ~ 2 / V2 ^ 2. 





X J dr] 

-1 

X 2F1 



■ ■ ■ in total {£ ~ l)/2 steps , 
until reaching 2^1(0, b; c; r/^) 



(73) 
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The coefficients ai{£, k) of the terms {£,Q—'if'Y^'^^^ are pro- 

/ a — I i 3 \ 

portional to 2^1 ( ^ h A;, - + 1 + fc; - + fc; 1 j and can 

be determined analytically. The terms formed by powers 
of ^0 can be summed up in closed form, and it turns out 
that they provide a very good approximation for the value 
of /^(eo): 



^^(^0 > 1) 



3^0 



(£-3)/2 

E 

s=0 



n 

.m=0 



2m + 3 



n 

.r?i.=0 





I 




(- 


i 








V m J 




- m j 




^ 2 " 






" 2 " 





2s+5 



2s + 5 



ai 



X 3F2 



3^0 



£2 + £ - 2 
45 



,3 £ ^ i 
l'2-2'^+2 



5 Zl..2 

2 ' 9 f ' s;o 



(74) 



Since Eq. ((7^ . which leads to a polynomial of order 2 — 1 

5 = £ + 2 [see the first equality in Eq. ([71|) ]. predicts that 
at large £ and ^0 > 1 the behavior is dominated by the 
term ^q"*"^, it obviously represents an approximation which 
breaks down for very large £ or for not small enough val- 
ues of ^0 ~ 1- However, this breakdown of the approxi- 
mation in Eq. (|74p can be easily monitored by checking 
that the result obeys the bound given by Eq. ([70)) . The 
approximation can be systematically improved, if needed, 
by including the terms ai(£, fc)(^g — 7^^)3/2+'^ (which can 
be computed analytically) as they become relevant. For 
?o - 1 < 10"'* and £ < 201 {£ odd), we did not have to 
add any such corrections to Eq. (fM]). 

The closed form of Eq. ((7^ significantly simplifies the 
numerical study of the behavior of Ie{£,o > 1) (with £ odd). 
It strongly supports (by extrapolating the results within 
the range where it provides an accurate approximation) 
that Jf^oo(Co) and that sign[/f(^o)] = sign(af) = 
(— 1)~2-. In combination with Eq. (p^ . this shows that 
approximating Ii{^o 1) by Eq. (j74p in the series repre- 
sentation of Eq. ((66)) leads to a general form of the term 
in the series for /(^o) which has alternating sign and van- 
ishes with increasing £, which ensures that as expected the 
series is convergent. 

We define the partial sum /l(^o), with L odd, of the 
series in Eq. ((66l) as 



\e=l, e odd ^^vsu; 



L+2 



+ E i^^+^)7^MCo)\ , (75) 

fcl, £ odd ^^(^0^ 

such that the oscillatory behavior induced by the alter- 
nating signs of successive terms is damped. Guided by the 
logarithmic divergence of Qi{£,a — > 1^), by the fact that 
the series defining the pref actor /(I) [Eqs. ([55)) and (|57)) ] 



turned out to be convergent (although vanishing) , and by 
the results for the phoretic velocity of a rod-like particle 
in Refs. [TTll^. we make the ansatz that /(^o) behaves as 



/(6>i) = (eo-i)in(eo-i)/(eo) 



(76) 



and aim at determining numerically the function /(^o)- 

In Fig. [5] we show the dependence of 
/l(Co) = fL{^o)/r{£.o) on L for various values e := ^0 - 1, 
IQ-M g iQ-i rpj^g convergence of the partial sums 
is clear, and this allows us to extract the values /(^o) as 
the corresponding constant plateau values of each of the 
curves /l(^o)- The numerically determined /(^o)j for clar- 
ity shown in the inset of Fig. [5] as a function of — log]^o(^)i 
highlights the following important features: (i) /(^o) < 
for all ^0 ^ 1, and (n) the data set {— logiQ(e), /(^o = 
1 + e)} is very well fitted by a power law [dashed line in the 
inset in Fig.[5] leading to /(^o 1^) — -0-7 [- logio(Co- 
l)]-0-9^-1.5 [-ln(eo-l)]-"-^. 



100 



200 



-0.2 



-0.4 



r 




-16 10 



-32 



10 



-64 







20 40 60 

» » Q Q S e ® ' 



-0.2 



■ l0g,„( E ) 



Fig. 5. Dependence of /z,(^o) (see the main text) on L 



(for L odd) for e 



Co 



10" 



10" 



10" 



and 



10" . The inset shows the resulting /(^o) (points) obtained 
from the extrapolated /z,->cx)(Co), as well as the power law fit 
— 0.7x [— logj^(,(e)]"'^ '' (dashed line), as a function of — logjQ(e). 
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